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ABSTRACT 

We present an analysis of the Far Ultraviolet Spectroscopic Explorer (FUSE) 
spectra of the little-known southern nova-like BB Doradus. The spectrum was 
obtained as part of our Cycle 8 FUSE survey of high declination nova-like stars. 

The FUSE spectrum of BB Dor, observed in a high state, is modeled with an 
accretion disk with a very low inclination (possibly lower than 10°). Assuming 
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was operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985 
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an average WD mass of O.8M leads to a mass accretion rate of lO _9 M /yr and 
a distance of the order of ~650 pc, consistent with the extremely low galactic 
reddening in its direction. The spectrum presents some broad and deep silicon 
and sulfur absorption lines, indicating that these elements are over-abundant by 
3 and 20 times solar, respectively 

Subject headings: accretion, accretion disks — novae, cataclysmic variables - 
stars: abundances — stars: individual (BB Doradus) - - stars: nova-likes - 
stars: white dwarfs 



1. Introduction 

1.1. Cataclysmic Variable Nova-Likes 

Cataclysmic variables (CVs) are close binaries in which the primary, a white dwarf 
(WD), accretes matter and angular momentum from the secondary, a main sequence star, 
filling its Roche-lobe. The matter is transferred by means of either an accretion disk around 
the WD, or an accretion column, when the WD has a strong (~ 10 MG) magnetic field. 
Ongoing accretion at a low rate (quiescence) is interrupted every few weeks to months by 
intense accretion (outburst) of days to weeks (a dwarf nova accretion event), and every few 
thousand years by a thermonuclear explosion (TNR; the classical nova event). CV systems 
are divided into sub-classes according to the duration, occurrence and amplitude of their 
outburst: dwarf novae (DN) are found mostly in the low state; while nova-likes (NL) exhibit 
the spectroscopic and photometric characteristics of novae between or after outbursts but 
have never had a recorded outburst. NLs form a less homogeneous class and are divided 
into sub- types depending on certain properties, such eclipsing systems, magnetic systems, 



or systems that go into unexpected low states (IWarner II1995I ). Non-magnetic NLs can be 
divided into UX UMa systems which remain in a high optical brightness state or "permanent 
outburst state", and the VY Scl systems (anti-dwarf novae) which experience unexpected 
low states when the optical brightness plummets. 

CVs are believed to evolve from long period (~2 day) toward short period (~1 hour) as 
the mass ratio q decreases, crossing the 2-3 hour 'period gap', where few system are found. 
Theory predicts that the period of those systems havin g a sufficiently low q wi ll increase as q 



continues to decrease, the so-called 'period bouncers' ([Patterson et al. 1 120051 ) - The driving 
mechanism behind the mass transfer (M) and hence, the evolution of CVs is believed to 
be angular momentum loss, dominated by magnetic stellar wind braking at periods above 3 
hours and by gravitational wave emission below this period. The secondaries are stripped 
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to < 0.08M Q on a timescale of only 1-4 Gyr after they form as close binaries above the gap. 
Consequently, in the lifetime of the galaxy, the vast majority of CVs should have evolved 
to a period minimum near 80 min and now hav e degenerate brown dw arf-like secondaries 



(IHowell et al. II2001I ). However, recent models (lAndronov et al. 1 120031 ) suggest (i) much 
lower angular momentum loss rates, so it takes 10-12 Gyr to even reach a 3 hr period, and 
(ii) multiple evolutionary tracks yielding different populations of CVs above and below the 
3 hour orbital period. 

The accretion disks and underlying accreting WDs in CVs can provide crucial clues and 
constraints on these evolutionary scenarios. It is important to derive M from observations, 
as the mass transfer rate are tied to the angular momentum loss. As to the WDs, they 
are the central engines of the observed outbursts, either as potential wells for the release of 
accretional energy or as the sites of explosive TNR shell burning. In addition, the differ- 
ences between the temperatures, rotation, and chemical abundances of CV WDs and single 
isolated WDs provide clue s as to the effects of accretion, diffusion, and long term heating 



and evolution (ISion 1119951 ). 

With a wavelength range from 905-1 187AA, FUSE covers that part of the far-ultraviolet 
(FUV) spectrum where the hotter CV component is dominant. For nova-likes, the dominant 
component in the FUSE range is usually the accretion disk. Additional possible components 
are the hot white dwarf, the boundary layer or accretion belt, and other hot regions on or 
close to the WD. Many CV FUV spectra are modeled with 2 components, usually a WD 
plus a disk. For each system, one usually derives from the FUSE spectrum the temperature 
of the WD (T e ff), the gravity, the rotational velocity (V ro tsmi), the chemical abundances, 
the mass accretion rate, the inclination, and the distance to the system. 

It has therefore been very important to try and obtain FUSE data for more CV systems 
in order to populate the T e jj-Period and the M- Period parameter spaces. There is a need to 
enlarge the sample to obtain an accurate global picture of the systems above and below the 
period gap. For that reason, we proposed to observe with FUSE of a set of 16 high declination 
dwarf novae (our Cycle 7 survey of DNe) and 16 high declination nova-likes (our Cycle 8 
survey of NLs - all chosen from the on-line CV catalog of Downes et al). Unfortunately, our 
FUSE NL survey was cut short due to the fatal failure of the reaction wheel of the telescope. 
Of the 16 targets, only two were observed: B B Dor and P831-57. The analysis of the FUSE 
spectrum of P831-57 is presented elsewhere ( Barrett et al. II2008I ). In this work we present 



the result of the FUSE spectroscopic analysis of BB Dor. 
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1.2. The Nova-like BB Doradus 



BB Dor (also known as EC 05287-5847, an object from the Edingburgh-Cape Survey 
(IChen et al. 1120011 )) was spectroscopically id entified as a CV in December 1987 and was seen 
fainter than V ~ 16.5 until November 1992 (IChen et al. Il200ll ). Since that time it seems to 
be in bright state ( V ~ 14.6 — 13.6). A first (extremely uncertain) estimate of its period by 
Chen et al. I (120011 ) put it right in th e middle of the p e riod gap with a period of 0.107 day or 
2.57 h. However, a 45 day coverage (IPatterson 1120021 ; iPatterson et al. 1 120051 ) corrected this 
value to 0.14923 d or 3.58 h, as expected for VY Scl stars which usually have perio ds between 
~ 3 and ~ 4 h. BB Dor exhibits superhumps with a period excess of e =9.4% (IPatterson 



20021 ) — it has the largest kno wn fractional superhum p of all systems and therefore provides 



a calibrating point at large e (IPatterson et al. II2005I ). BB Dor is a VY Scl nova-like and so 



far no magnetic act i vity has been found . The system also exhibits quasi-periodic oscillations 
( IChen et al. Il200ll ). I Chen et al. I (120011 ) also noted that due to its very narrow Balmer lines 
(with wings extending ~ 650km s _1 from the line center) BB Dor is most likely a low- 
inclination system. From the Austral Variable Star Observer Network (AVSON), it appears 
that BB Dor has been around V ~ 13.6 — 14.6 mag (with an error < 0.1 mag) for the last 
couple of years and can vary as much as 0.4 mag within one hour. The reddening value E(B- 
V) towards BB Dor is not known. However, the Galacti c reddening value in t he direction of 
BB Dor as inferred from the 100/im dust emission map (ISchlegel et al. Ill998l ) is very small: 
E(B-V)=~0.03. Since the dust map gives the maximum dust extinction, we deduce that BB 
Dor must have a negligible reddening (E(B — V) ~ 0.01 at most) which should not affect 
the results; consequently we do not deredden its FUSE spectrum. 



2. Observations 

2.1. The FUSE Spectrum of BB Dor 

BB Doradus was observed with FUSE on 09 July 2007 at 02h51m09s (UT), and 3 
exposures were obtained (H9030301 - exposures 001, 002, 003). The FUSE spectrum was 
obtained through the low resolution (LWRS; 30"x30") aperture in TI ME TAG mode. Th e 



data were processed with the last and final version of CalFUSE v3.2.0 (IDixon et al. 1 120071 ). 



A look at the count rate graph (generated by CALFUSE) indicates that the data for each 
exposure were collected mainly during the first 1000 s, followed by a jitter lasting about 1,500 
s. Toward the end of each exposure, the count rate dropped to zero. Therefore, while the 
co-added exposures for the 8 individual segments after processing by CalFUSE give a good 
exposure time (after screening) of between about 8750 s (segments 1) and 9750 s (segments 
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2), the count rate plot indicates that the total good collection time ranges between 3350 s 
(segments 1) and 3650 s (segments 2). We therefore weighted these segments accordingly, 
multiplying the flux (of the order of 1 x 10~ 13 ergs s _1 cm _2 A -1 in the on-line preview) by a 
factor of 2.61 (segments 1) and 2.67 (segments 2). After that we processed the data as usual 
(see below). 

The spectral regions covered by the spectral channels overlap, and these overlap regions 
are then used to scale the spectra in the SiCl, SiC2, and LiF2 channels to the flux in the LiFl 
channel. The low sensitivity portions (usually the edges) of each channel are discarded. In 
the present case the SiC2b spectral segment was very noisy and we discarded it too. We also 
carried out a visual inspection of the FUSE channels to locate the worm and we manually 
discarded these portions of the spectrum affected by the worm. We combine the individual 
exposures and channels to create a time-averaged spectrum weighting the flux in each output 
datum by the exposure time and sensitivity of the input exposure and channel of origin. The 

final product is a spectrum that covers the full FUSE wavelength range 905 1187AA. 

Because we disregarded the right edge of both the SiCla and SiC2b segments, there is a gap 
between the LiFla and LiF2b segments (~ 1082 - ~ 1087AA). 

The FUSE spectrum of BB Dor is shown in Figure 1. Note that the FUSE exposures 
of BB Dor consist entirely of observations carried out during the DAY. The small NIGHT 
exposure was ignored because of jitters and other problems. We suspect that most of the 
emission lines are due to helio- and geo-coronal emissions. 

2.2. The FUSE Lines 

The FUSE spectra of CVs exhibit mainly broad absorption lines from the accretion disk 
and WD, as well as sharp absorption lines from circumstellar (or circumbinary) material 
and the ISM. Emission lines from the hotter regions are usually broadened, depending on 
the inclination angle i, due to the large Keplerian velocity in the inner disk. Sharp emission 
lines from air glow (geo- and helio-coronal in origin) are also present, due mainly to the 
reverberation of sunlight inside the FUSE telescope during DAY time observations. The 
present FUSE spectrum of BB Dor presents such a complexity, and all the lines are listed in 
Table 1. 

Absorption Lines. The main characteristic of the FUSE spectrum of BB Dor is the broad 
Ly(3 absorption feature due to either the exposed WD, the disk at low inclination angle 
i, or possibly both (see the Results section). Absorption features from the higher orders 
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of the Lyman series are also clearly visible, indicating a temperature T~ 35, 000K. The 
other main absorption features expected at this temperature are the Cm (1175 A), Sim 
(wll08-1114AA and «1 140-1 144AA), Si IV («1120-1130AA) and the Cn (lOlOA). The 
spectrum, is however, also marked by some high-order ionization absorption lines such as 
Siv (1063 & 1073A), Si IV (1066.6A), and the O vi doublet. The Siv (1006A), Cn (lOlOA), 
and Siv (1099 A) absorption multiplets are all unresolved and shallow. However, we have 
indicated them on the figures, as they appear deeper in the modeling. On the other side the 
Siv (1063A & 1073A) and Sim (~ 1010A) absorption lines appear deeper in the observed 
spectrum than in the modeling. The presence of the high-order ionization species (such as 
the oxygen doublet) indicates the presence of a hotter absorbing component above the main 
FUV emitting region. The Pll line (~96lA) has been marked in the spectrum, though it is 
most probably a feature due to the noise. All the the lines are listed in Table 1 with their 
wavelengths. Most of the broad absorption lines and features associated with the source are 
redshifted by about 0.6-0.9A, corresponding to a receding velocity of 200 ± 40km s^ 1 . It 
seems very likely that this redshif is due to the orbital motion of the WD, however, we do 
not have enough time tag data to verify this assumption. 



Emission Lines. The spectrum of BB Dor exhibits sharp emission lines including all the 
orders of the Lyman series, Niv (~923A), S vi (933.5A & 944.5 A), Cm (977A), Hell (992A), 
and He I (1068 A). Some of these lines seems to be due to helio-coronal emission (sunlight 
reflected inside the telescope) that contaminates the SiC channels, in particular the H I lines, 
the He I & Hell lines, the Cm (977A) line, and the Ovi doublet. The segments of the 
spectrum from the LiF channels overlapping the SiC Channels down to about 1000A do not 
show any oxygen lines, and we therefore did not include these spectral regions of the SiC 
channels. However, other emission lines may be solely from the source. In many NLs, the 
emission lines from the dis k are broadened by the Keplerian velocity and are consequently 



easily identified. However, IChen et al. I (120011 ) have detected very narrow Balmer lines and 



concluded that BB Dor is most likely a low-inclination system. Because of this, the identifi- 
cation of the emission lines from the source is not trivial. This issue is further complicated by 
the lack of NIGHT exposure and the very low signal-to-noise ratio, which makes it difficult 
to analyze the two-dimensional image of the spectrum for the effect of scattered sunlight on 
the emission lines. 

Since the system is a NL in a high state, we expect at first the Niv (?^923A), Svi 
(933.5 & 944.5 A), and Cm (977A) emission lines to be from the source with some possible 
helio-coronal contamination mostly affecting the carbon emission lines. The Niv (923.06A) 
emission line is contaminated with the Hi (923. 15A) line; the Niv (922. 52A + 924.28A) 
lines are not. However, the complete absence of narrow emission lines from the oxygen 
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doublet and the Cm (1175) seems to indicate otherwise — i.e. all the sharp emission lines 
might be helio-coronal in origin. To further check this possibility we measured the narrow 
emission line relativ e intensities below lOOOAand compared with those of the solar spectrum 
( ICurdt et al. Il200ll ). We found a good agreement with the solar-disk average quiet Sun data, 
except for the SVI 944.5 which has a higher relative intensity in the spectrum of BB Dor. 
This could be due sun spot activity during the FUSE observation of BB Dor, as this line 
intensity increases by a factor of 30 inside sun spots while the other lines increase at most 
by a factor 3. We therefore conclude that all the sharp emission lines in the FUSE spectrum 
of BB Dor are due to sun light reflected inside the telescope. 



3. Spectral Modeling 

We created a grid of models of synthetic spectra of WDs and accretion disks for different 
values of the WD temperature T e ff, gravity Log(g), projected rotational velocity V rot sin(z), 
inclination i, mass accretion rate M, and abundances. We then ran a y 2 fitting program to 
find the best fit for (i) a single WD component, (ii) a single accretion disk component, and 
(hi) a combined WD + accretion disk. We describe how we generate these synthetic spectra 
and how we perform the fitting. 

We create the syntheti c model spectra for high-gravity stellar atmospheres using codes 
TLUSTY and SYNSPEC0 Jrlubenv Ill988l : lliubenv fc LanTlll995h . Atmospheric structure 



is computed (using TLUSTY) assuming a H-He LTE atmosphere; the other species are then 
added in the spectrum synthesis stage using SYNSPEC. We generate photospheric models 
with effective temperatures ranging from ~20, 000-50, 000 K in increments of 1,000 K. We 
chose values of Log(g) ranging between 7.5 and 9.0. We also varied the stellar rotational ve- 
locity V ro tsin(i) from 100km s _1 to 1000km s -1 in steps of 100km s™ 1 (or smaller if needed). 
The WD rotation (V rot sin(ij) rate is determined by fitting the WD model to the spectrum 
while paying careful attention to the line profiles in the FUSE spectrum. We do not carry out 
separate fits to individual lines but rather try to fit the lines and continuum simultaneously 
while paying careful attention to the absorption lines. It is important to note that the depth 
of the absorption features depends not only on the abundances but also on the rotational ve- 
locity. Increasing the rotational velocity reduces the central depths of the absoption features, 
thus reducing the abundance. However, the widths of the absorption features also increase 
with increasing rotational velocity. As a consequence, abundances and rotational velocity 
effects are interwined and cannot always be easily separated without additional information 



i 



http: / /nova.astro.umd.edu] TLUSTY version 200, SYNSPEC version 48 
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on either parameter. For any WD mass there is a correspon ding radius, or equiva l ently one 
sin gle value of Lo q (q) (e.g. see t h e mas s radius relation from lHamada fc Salpeter I (Il96ll ). or 



sec 



Wood I (11990 ): IPanei et al 



(120001 ) for different composition and non-zero temperature 
WDs). Therefore, by scaling the theoretical spectrum to the observed one, we obtain the 
distance to the system. 

We model accretion disk spectra by first assuming that the disk is made of a collection 
of annuli, where each annulus has a temperature T(r) and gray ity Log(g(r)) given by the 
standard disk model (IShakura fc Sunyaev I Il973l ; iPringle I Il98ll ). for a given central mass 
Mwd and accretion rate M. A variant of the code TLUSTY is then used, TLDISK, that 
generates an atmosphere model for each annuli that is then used as input for the code 
SYNSPEC. The contribution of all the annuli are then combined using the code DISKSYN, 
and a final spectrum is o btained for any given inc lination angle. A detailed explanation of 
the procedure is given in I Wade fc Hubeny I (119981 ). I n the present work we do not use the 
grid of synthetic accretion disk spectra tabulated by IWade fc Hubeny I (119981 ). instead we 
generate them. This allows us to compute disk spectra assuming no n-solar abundances and 
for any inclination angle (the disk spectra of I Wade fc Hubeny I (119981 ) have solar abundances 
and have been generated for given value of the inclination i). 

Before carrying out a synthetic spectral fit of the spectra, we masked portions of the 
spectra with strong emission lines, strong ISM absorption lines, detector noise, and air glow. 
The regions excluded from the fit are in blue in Figures 1 and 2. 

After having generated grids of models for the FUSE spectrum of BB Dor, we use FIT 
(IPress et al.lll992l ). a x 2 minimization routine, to compute the reduced xl {x 2 P er number of 
degrees of freedom) and scale factor (which gives the distance) for each model fit. While we 
use a x 2 minimization technique, we do not blindly select the least x 2 models, but we also 
examine the models that best fit some of the features such as absorption lines and, when 
possible, the slope of the wings of the broad Lyman absorption features. 

Initially we generate solar abundances models, and when a good fit is found, we start 
varying the chemical abundances of C, N, S and Si, to fit the absorption features of the 
spectrum. In particular, the carbon abundance was set using the Cn (1010A) and Cm 
(1175 A) multiplets; the sulfur abundance was set using the S IV (1063A & 1073A) lines; and 
the silicon was set using the Si IV (1067A , 1023A , 1028A) and Sim (~1110A , «1138-1146A) 
lines. 
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4. Results 

The data obtained by the AVSON implies that BB Dor has been in a high state around 
V ~ 13.6 — 14.6 for the last couple of years. We therefore expect the FUSE spectrum to 
be dominated by flux from the accretion disk with a relatively high M. However, in our 
modeling we follow a systematic approach which consists in fitting the following: (i) a single 
WD; (ii) a single accretion disk; and (iii) a WD + accretion disk composite. 

4.1. White Dwarf. 

Since we do not have any information on the mass of the WD, we look for all the best-fit 
models in the Log(g) versus T e ff plane. Namely, for each assumed value of Log(g) we vary 
the temperature to find the best-fit model. As expected we find that the temperature is 
somewhere between 34,000K (for Log(g) = 7.5) and 40,000K (for Log(g) = 9.0), with a 
distance between 377 pc and 135 pc, respectively (see Table 1). The least xl is obtained for 
Log(g) = 8.0 — 8.65. We chose the intermediate value Log(g) = 8.3 model, with T=37,000K 
and d=247 pc and solar abundances, to illustrate our results in Figure 1. We then further 
improve the fit by varying the abundances. We find that in order to better fit the sulfur 
and silicon lines, we have to set the sulfur to 20 times its solar abundance, and the silicon 
to 3 times, while all the other species are kept at solar abundances. The xl value decreases 
from 0.3348 to 0.3129 (by about 7%). While some of the lines are better fitted (such as Siv 
1063A & 1073A; Sim ~1010A & SilV -1025A), the Siv (1006A & HOOA) are far too deep 
(see Figure 2). 

4.2. Accretion Disk. 

Next, we fit the solar abundances disk models. We find that a low inclination is needed 
in order to match the absorption lines, and in our models we initially set i=5, 8, 12, and 18 
degrees. Again, because we have no information on the WD mass, for each value of Log(g) 
(ranging between 7.5 and 9.0) we vary the mass accretion rate between 10~ 10 ' 5 M o yr~ 1 and 
lO _8 -°M yr" 1 to find the best-fit models. Again we find that the best-fit models are around 
Log(g) = 8.3, but the difference between the xl values is not significant, and the improvement 
over the single WD models is also only marginal (of the order of ~1% in the xl value). All 
these models are presented in Table 2. Next, we decide to improve the disk model by varying 
the abundances, though we do not expect the fit to improve much because of the Keplerian 
velocity broadening of the disk. However, we find that setting the silicon and sulfur to 3 
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and 20 times solar, resp., actually reduced the xl value by 10%, more than for the single 
WD model. This is because the the Siv (1006A & llOOA) lines have better fits (Figure 
3). This best-fit single disk model has % = 8°, M = lO _9 M yr _1 , a distance of 665 pc and 
Xl = 0.3000. 



4.3. Composite Model: WD + Accretion Disk. 

Last, we fit composite WD+Disk models. Since the number of models increases expo- 
nentially when we change Log(g), T e ff and M, we restrict our search using our best fit for 
Log(g), namely 8.3. We also set the sulfur to 20 times solar and silicon 3 times solar. This 
procedure reduces the number of free parameters to 3 (M,T e ff,i). 



Low Inclination. Since the system is believed to have a low inclination angle (IChen et al. 



20011 ). we generate low-inclination models (i = 5°, 8°, 12°, & 18°). Not surprisingly, the 
best fit model has % = 8° and M = 10 -9 M Q yr -1 , but this time the WD has a temperature of 
32,000K and provides only 10% of the flux, while the disk provides the remaining 90%. Such 
a model again brings an insignificant improvement in the value of xl an d it is clear (Figure 
4) that it is barely distinguishable from the best-fit single disk model. From the point of view 
of the physics, this model is preferred because, if the system has a low-inclination angle, then 
the emission from the WD must contribute to the flux. Models with a WD temperature T < 
30, 000K had a WD contribution of only a few percent of the total FUV flux, and could not 
be distinguished from the single disk models. If the WD has a temperature T e ff <30,000K 
it will not be detected while the system is in a high state (with M = lO^M^r" 1 ). These 
results implies that the contribution from the WD is not very large and that the temperature 
of the WD must be < 32, 000 K. 



High Inclination. From the sharp Balmer emission lines, I Chen et al. I (120011 ) suggest the 
system has a low inclination. However, we cannot confirm the inclination directly from the 
emission lines of the FUSE spectrum, as all the sharp emission lines in the FUSE spectrum 
are of helio-coronal origin. Our low-inclination angle single disk models provide a slightly 
better fit than the single WD models (and a much better fit than the high-inclination single 
disk models), however, for completeness we include here the results from the composite WD 
+ disk model fits when the assumption on the inclination angle is relaxed. 

We search for the best-fit WD+disk models (assuming Log(g) = 8.3) in the parameter 
space T e ff versus M including all inclination angles, and find that the models with an 
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intermediate inclination (i = 18°, 41°, 60°, & 75°) do not provide the best fit. The best-fit 
model has i = 80°, and reflects a situation in which the WD is dominant with T = 37, 000K 
and contributes 3/4 of the total flux, while the disk has M = lO _9 M yr _1 and contributes 
only 1/4 of the flux. The distance obtained from this model is 246 pc and xl — 0.2873. This 
is the least x 2 value we obtained from all our models. This model is presented in figure 5. 
Similar results were obtained when assuming Log(g) = 8.65, but with a slightly lower mass 
accretion rate (Table 2). 



5. Discussion and Conclusion 

BB Dor is a little-known southern NL, and, consequently, both the distance and the 
mass of the WD are unknown, which implies a larger uncertainty in the results. In addition 
the FUSE spectrum is definitely of poor quality. In theory, a fine tuning of the temperature 
(say to an accuracy of about ±50 K) and mass accretion rate can be carried out by fitting 
the flux levels such that the distance to the system (when known) is obtained accurately. 
However, the fitting to the distance depends strongly on the radius (and therefore the mass) 
of the WD. We discuss some additional restrictions that we use to constrain the properties 
of the system. 

On the basis of the least x 2 v i ^ e best model is the high inclination WD+disk composite. 
However, (see Figure 5), because of the high inclination, the disk contributes a rather flat 
component all the way into the shorter wavelengths. At such a high inclination, one would 
not expect the WD to dominate the flux, but rather the WD would be almost completely 
masked by the swollen disk; actually the system would likely be observed to undergo eclipses 
- but none have been observed. Also, since BB Dor can vary as much 0.4 mag in one hr, this 
means that the light cannot be dominated by the WD, and the disk must be contributing 
at least 40% of the light. In other words, while the least x 2 indicates that the best fit 
is a WD with a rather flat disk component, there are other indications that this cannot 
be correct. Actually, the WD+disk high inclination model strikingly resembles the second 
com ponent observed in th e FUSE spectra of some dwarf novae during quiescence (e.g. VW 



Hyi, iGodon et al. I (120041 ) ). This is likely a caveat in the state-of-the-art spectral modeling, 
rather than an indication of a physical link between the the spectrum of BB Dor and that 
of a DN in quiescence. The need for improved modeling also stems from the difficulty in 
producing a model that fits low and high-order ionization lines at the same time. For BB 
Dor it is possible that the S IV 1063A & 1073A absorption lines form in the same hotter 
region/layer where the Ovi doublet form; while the Cm and Sim lines form in a cooler 
region/layer where the Lyman series (and continuum) form. This is similar to the HST/STIS 
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spectrum of TT Crt (jSion et al. 1120081 ) which exhibits a rich variety of absorption lines from 
different ionization stages, suggesting line formation in (at least) two different temperature 
regions. 

In order to reduce the size of the domain for which we have best-fit models in the 
parameter space, we use the infrared magnitudes J, H, & K from the Two Mic r on All Sky 
Survey (2MASS) to assess the distance to BB Dor as prescribed by IKnigge I (120061 ) for 
systems with P < 6 h. The IR data were collected on 9 November 1999, at a time when 
BB Dor was in a high state with a visual red magnitude R=14.60 and a blue magnitude 
B=13.90 (whereas B~16.5 in the low state). The 2MASS IR apparent magnitudes are 
J=14.322, H=14.089, and K=14.053. For a primary star M wd = 0.7 5M m and period o f 3.559 



O.25M ([Patterson et al 



h (corresponding to BB Dor), the donor star mass is M 2 
suggest M 2 = 0.256M Q ) and the IR absolute ma gnitude e s timat es are Mj 
and Mk = 6.63. Inserting these into eq.(15) of IKnigge I (120061 ) gives a distance of 235, 274 



(120051 ) 
7.47, M H = 6.90 



& 305 pc, respectively. These distances are typically underestimated by factors of 2.05, 1.86 
and 1.75 for the J, H and & K bands giving distances of 482, 510, and 534 pc, respectively , 
assuming the donor star contributes only ~ 1/4 — 1/3 of the total IR flux (IKnigge II2006I ). 
The distance to BB Dor is therefore certainly larger than 300 pc, and most likely in the 
range of ~ 500 pc. Since BB Dor was observed in a high state (with 2MASS), it is likely 
that the infrared flux of the donor star is on the lower side, namely even less than 1/4, in 
which case the 500 pc itself is only a lower limit. Using this restriction on the distance, we 
see that we can safely reject the single WD models (with d < 400pc), the high inclination 
models (with d < 250pc), and the large mass Log(g) = 9.0 models (with d = 361pc and 
less) . 

On the basis of the distance alone, the low inclination disk and WD+disk models are 
the best fits. These are also the best fit on the basis of the optical spectrum (low inclination) 
and visual magnitude (high state/disk). 

We are now able to summarize the basic characteristics of the system and its WD. 
We are confident in the fact that the system must have a WD with Log(g) = 8.3 ± 0.3, 
T e fj = 32,000K or lower (disk+WD composite models), a mass accretion rate of the order 
of M — 10 _9 M Q yr _1 (disk and disk+WD models) and a distance somewhere between 500 
and 700 pc (choosing the least \ 2 models agreeing with the distance). From the best-fit 
single disk models alone we find that the inclination angle must be very small (~ 8°) to fit 
the absorption lines. There is an indication that the sulfur and silicon are over abundant 
and it is also possible that some of the high-order ionization lines form in a layer in front of 
the disk and the WD. 



Although the FUSE spectrum of BB Dor is rather poor and multiwavelength data are 
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limited, we have narrowed down the region in the parameter space (M, P, T e ff) for this 
system. BB Dor is now the sixth VY Scl NL variable with a temperature estimate for 
its WD. With a temperature of < 32, 000 K, the WD of BB Dor marks the lower end 
of the temperature distribution for WDs of VY Scl NL variables, which, so far, for the 5 
other systems, ranged between 40,000K and 47,000K jHamilton fc Sion Il2008h . The FUSE 
spectrum of BB Dor clearly shows the drop in flux in the shorter wavelengths, the signature of 
a moderate tempera ture when compared to the FUSE spectra of these other VY Scl systems 
(see e.g. V794 Aql, (IGodon et al. 1120071 )). Even if we consider the single WD best-fit model, 
the WD temperature of BB Dor only reaches 37,000K (and the 40,000K Log(g) = 9.0 model 
gives an unrealistic close distance of 135pc in disagreement with the expected IR emission 
from the secondary). In that respect, the data point in the (T e fj, P) plane for BB Dor falls 
under the region of the VY Scl N L variables, seemi ngly into the region of the parameter 
space populated with dwarf novae (jSion et al. II2008I ). However, as shown in Figure 6, the 
data point for BB Dor does not seem to stand apart, and the separation between VY Scl 
systems and DN can be made easily by drawing a diagonal line. In that respect, it is actually 
V794 Aql that divides the two regions of the graph as indicated by the slanted line. 



We thank the anonymous referee for a prompt, concise and pertinent report which helped 
improve the analysis of the lines. Special thanks to the Austral Variable Star Observer 
Network (AVSON), and in particular to Berto Monard for providing us with the visual 
magnitude of BB Dor between 02 Sep., 2005 and 30 Dec, 2007, and to Peter Nelson for a 
quick look V band photometry at BB Dor on 28 Jan., 2008. PG wishes to thank Mario Livio, 
for his kind hospitality at the Space Telescope Science Institute, where part of this work was 
done. The infra-red J, H, & K magnitudes of BB Dor were retrieved from the online archival 
data of the 2MASS at IPAC Part of this publication makes use of the data products from the 
Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the National Science Foundation. 
This research was based on observations made with the NASA-CNES-CSA Far Ultraviolet 
Spectroscopic Explorer. FUSE is operated for NASA by the Johns Hopkins University 
under NASA contract NAS5-32985. This research was supported by NASA FUSE grant 
NNX07AU50G to Villanova University (P. Godon). This work was also supported in part 
by NSF grant AST05-07514 to Villanova (E.M. Sion). 
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Table 1. FUSE Line Identification 



Ion 


Kest (,AJ 


Kbs (A) 


Comments 61 


H T 

11 1 


918 13 

JlU. io 


918 




c/e 




Q1 Q 3^ 


919 3 




c/e 




Q2D 96 


920 9 




c/e 


Ntv 

IN IV 


999 ^9 


922 3 




s/e 




993 06 


923 




s,c/e 


H T 

11 1 


993 1 ^ 


923 2 




c/e 


N TV 

IN 1 V 




Q9S 9 




c/e 




994 28 


924 2 




c/e 


H T 

11 1 


996 93 


926 2 




c/e 






qso 7 




c/e 


S VT 


933 ^0 


933 4 

kJKJKJ .TT 




c/e 


H T 

11 ± 


937 8D 

<_/*J 1 .Uu 


937 8 




c/e 


S VT 


Q44 ^0 

• 7 tit: .UU 


944 4 




c/e 


H T 

11 1 


Q4Q 74 

C/TrC/ . I TT 


Q4Q 7 




c/e 


Nt 

1 N 1 


9^9 4D 


92^ ^ 




ism/a 




49 




















Q"i3 Q7 












9^4 9 




ism/ a 


P TT 

J. 11 


Qfil (14- 

Jul .Ut 


Qfi1 9 

Jul .Zj 




ism/a? 


H T 

11 1 


Q72 ^4 


979 5 




c/e 


Cm 


977.02 


977.0 




c/e 


Sim 


993.52 


993.2 




s?/a 




994.19 


994.9 




s?/a 




997.39 


997.1 




s?/a 


Siv 


1006.07 




s/a 


- unresolved 




1006.39 




s/a 


- unresolved 


Cii 


1009.85 




s/a 


- unresolved 




1010.08 




s/a 


- unresolved 




1010.37 




s/a 


- unresolved 


Hi 


1025.77 


1025.7 




c/e 
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Table 1 — Continued 



Ion 


Arest (A) 


\obs (A) 


Comments" 


Ovi 


1031.91 


1032.9 




s/a 


Cn 


1036.34 


1036.4 




ism/ a 


Ovi 


1037.61 


1038.3 




s/a 


Ol 


1039.10 


1039.3 




/ 

ism/a 


An 


1048.20 


1048.3 




/ 

ism/ a 


Siv 


1062.65 


1063.4 




/ 

s/a 


Si iv 


1066.60 


1067.1 




/ 

s/a 


An 


1066.66 


1066.8 




/ 

ism/a 


Siv 


1072.97 


1073.9 




s/a 




1073.52 


1073.9 




/ 

s/a 


Siv 


1098.36 


— 


/ 

s/a 


- unresolved 




1098.93 


— 


/ 

s/a 


- unresolved 




1099.48 


— 


/ 

s/a 


- unresolved 




1100.53 


— 


s/a 


- unresolved 


Sim 


1108.36 


1108.8 




s/a 




1109.94 


1110.6 




s/a 




~1113.20 


1113.8 




/ 

s/a 


Si iv 


1122.49 


1123.0 




/ 

s/a 




~1128.33 


1128.9 




s/a 


Ni 


1134.16 


1134.3 




c/a 




1 1 OA AO 






c/a 




1134.98 


1135.2 




c/a 


Si 


1145.10 


1145.0 




ism/a 


Hen 


1168.61 


1168.6 




c/e 


Cm 


1174.90 


1176.4 


s/a 


- unresolved 




1175.26 


1176.4 


s/a 


- unresolved 




1175.60 


1176.4 


s/a 


- unresolved 




1175.71 


1176.4 


s/a 


- unresolved 




1176.00 


1176.4 


s/a 


- unresolved 




1176.40 


1176.4 


s/a 


- unresolved 



-18- 



a The following abbreviations have been used: 
a - for absorption; e - for emission; c - for con- 
tamination (e.g. air glow; helio-coronal emis- 
sion); s - for source; and ism - for interstellar 
medium. 
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Table 2. Synthetic Spectra 



T 1 \ 

Log(g) 


wd 


V rot sin% 


[Si] 


[S] 


[Z] 


i 


Log(M) 


WD/ Disk 


d 


2 


Fig. 


(cgs) 


(10 K) 


(km/sj 


Solar 


Solar 


Solar 


(deg) 


(M /yr) 


(%) 


(P c ) 






7.50 


34.0 


400 


1.0 


1.0 


1.0 






100/0 


377 


0.3427 




8.00 


36.0 


400 


1.0 


1.0 


1.0 


— 


— 


100/0 


266 


0.3361 




8.30 


37.0 


400 


1.0 


1.0 


1.0 






100/0 


247 


0.3348 


1 


8.30 


37.0 


400 


3.0 


20. 


1.0 






100/0 


211 


0.3129 


2 


8.50 


38.0 


400 


1.0 


1.0 


1.0 






100/0 


193 


0.3345 




8.50 


38.0 


400 


3.0 


20. 


1.0 






100/0 


184 


0.3100 




8.65 


38.0 


400 


1.0 


1.0 


1.0 






100/0 


172 


0.3360 




9.00 


40.0 


400 


1.0 


1.0 


1.0 






100/0 


135 


0.3390 




7.50 






1.0 


1.0 


1.0 


12 


-8.0 


0/100 


878 


0.3346 




7.88 






1.0 


1.0 


1.0 


12 


-8.5 


0/100 


878 


0.3346 




8.30 






1.0 


1.0 


1.0 


8 


-9.0 


0/100 


693 


0.3297 




8.30 






3.0 


20. 


1.0 


8 


-9.0 


0/100 


665 


0.3000 


3 


8.65 






1.0 


1.0 


1.0 


5 


-9.5 


0/100 


496 


0.3332 




9.00 






1.0 


1.0 


1.0 


5 


-10.0 


0/100 


361 


0.3376 




8.30 


32.0 


400 


3.0 


20. 


1.0 


08 


-9.0 


10/90 


700 


0.2989 


4 


8.30 


37.0 


400 


3.0 


20. 


1.0 


80 


-9.0 


74/26 


246 


0.2873 


5 


8.65 


38.0 


400 


3.0 


20. 


1.0 


80 


-9.5 


75/25 


184 


0.2880 
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Fig. 1. — Solar abundance single WD model. The best-fit WD model assuming solar 
abundances (solid black line) is shown together with the FUSE spectrum of BB Dor (light 
grey/red line). The segments of the spectrum that have been masked before the fitting are 
in dark grey/blue. The WD has a temperature T = 37, 000K, a projected rotational velocity 
V rot smi = 400km/s, and Log(g) = 8.3. The distance is d =217 pc, and xl — 0.3348. 




Fig. 2. — Non-solar abundance single WD model. The best-fit WD model to the FUSE 
spectrum of BB Dor is shown where the abundances of sulfur and silicon have been set 
to: [S]=20 solar and [Si] =3 solar. The abundances of all the other elements have been 
set to solar. The WD has a temperature T = 37, 000K, a projected rotational velocity 
V rot smi = 400km/s, Log(g) = 8.3. The distance obtained is d =211 pc, and xt = 0.3129. 
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Fig. 3. — Non-solar abundance single disk model. The FUSE spectrum of BB Dor is shown 
together with a one of the best-fit synthetic disk models. The model has M = O.8OM , 
M = 10~ 9 M Q /yr, i = 8°, and gives a distance of d = 665 pc and a xl — 0.3000. The sulfur 
and silicon abundances have been set to: [S]=20 solar and [Si] =3 solar. The abundances of 
all the other elements have been set to solar. The low inclination is needed in order to fit 
the silicon, sulfur and carbon absorption lines. 
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Fig. 4. — Low-inclination non-solar abundance WD+disk model. The FUSE spectrum of BB 
Dor is shown together with a one of the best-fit WD+disk models (solid black line), assuming 
a low inclination angle. The model has M = 0.80M Q , M = 10 _9 M Q /yr, % = 8°, and gives a 
distance of d = 700 pc and a. xl = 0.2889. The sulfur and silicon abundances have been set 
to: [S]=20 solar and [Si] =3 solar. The abundances of all the other elements have been set 
to solar. The WD (dotted line) contributes 10% of the flux while the high-inclination disk 
(dashed line) contributes the remaining 90%. 
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Fig. 5. — High-inclination non-solar abundance WD+disk model. The FUSE spectrum of 
BB Dor is shown together with a one of the best-fit WD+disk models (solid black line). 
The model has M = O.8OM , M = lO _9 M /yr, and gives a distance of d = 246 pc and a 
xl = 0.2873. The inclination angle of the system is kept as a free parameter and the best-fit 
model gives % = 80°. The sulfur and silicon abundances have been set to: [S]=20 solar and 
[Si] =3 solar. The abundances of all the other elements have been set to solar. The fit of 
the silicon, sulfur and carbon absorption lines is obtained by having the WD component 
dominating the spectrum. The WD (dotted line) contributes 74% of the flux while the 
high-inclination disk (dashed line) contributes the remaining 26%. 
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Fig. 6. — Dwarf nova systems and VY Scl nova-like variables above the Period-gap. The 
systems have been plotted in the (T P ff,P) parameter space on a logarithmic scale. The data 
for the nova-likes are taken from lHamilton fc Sion I (120081 ) and references therein; the data 
for the dwarf nova systems are taken Ision et al.~ ( 2008h and the refere nces therein. For RU 
Peg we have used the more recent estimate from (IGodon et al. 1 120081 ) . The data point for 
V794 Aql is an estim ate of the WD temp erature based on an IUE spectrum while the system 
was in a lower state (IGodon et al. 1120071 ). While the data point of BB Dor is lower than the 
other VY Scl system there is still a clear separation between the VY Scl systems and the 
DN systems as indicated by the slanted line. 



